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Abstract. We present the results of the first N-body simulations of the Growing Neu- 
trino scenario, as recently discussed in |Baldi et af] ( |201 1[ ). Our results have shown 
for the first time how neutrino lumps forming in the context of Growing Neutrino cos- 
mologies are expected to pulsate as a consequence of the rapid oscillations of the dark 
energy scalar field. We have also computed for the first time a realistic statistical dis- 
tribution of neutrino halos and determined their impact on the underlying Cold Dark 
Matter structures. 



1. Introduction: the Growing Neutrino scenario 



The Growing Neutrino scenario has been recently proposed by Amendola et al. (2008) 
as a possible solution to the so called "cosmic coincidence" problem, i.e. the fact that 
the two main components that constitute our Universe, Dark Energy (DE) and Cold 
Dark Matter (CDM), happen to have comparable densities only around the present cos- 
mological epoch. The Growing Neutrino scenario is based on a direct interaction be- 
tween a DE scalar field <p, which is considered to be responsible for the late-time accel- 
eration of the expansion of the universe instead of the standard cosmological constant, 
and massive neutrinos. This interaction determines an exchange of energy between the 
two components such that the individual energy-momentum tensors of DE and massive 
neutrinos are not separately conserved, while only their sum obeys the standard conser- 
vation equation provided by General Covariance. This exchange of energy is encoded 
in the dependence of the average neutrino mass on the global evolution of the DE scalar 
field <p, according to the equation: 

m v = m v (t )e-^ , (1) 



where m v (?o) is a constant that represents the value of the neutrino mass at the present 
time and the parameter /? sets the strength of the coupling. The dynamics of the DE 
scalar field <p will therefore determine a time dependence of the neutrino mass, and 
for a given set of cosmological parameters the background evolution of the system 
can be computed. In particular, a specific feature of viable Growing Neutrino models 
(characterized by a large and negative value of the coupling /3, which is taken to be 
(3 = -52 for the realization considered in the present work) is the sudden growth of 
the neutrino mass at the time when neutrinos become non-relativistic particles. In fact, 
while for relativistic particles the effect of the coupling is suppressed by pressure terms 
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in the neutrino continuity equation, as soon as the neutrino equation of state parameter 
w v = Pv/Pv drops from its relativistic value of 1/3 to zero the coupling becomes active 
and introduces a steep minimum in the total effective scalar potential where the scalar 
field <p rapidly oscillates and eventually stops. 

The trapping of the DE scalar field in a steep local minimum of its effective poten- 
tial determines a sudden drop of the DE equation of state w$ towards the value of - 1 , 
and the DE scalar field therefore starts to behave as a cosmological constant. The tran- 
sition from relativistic to non-relativistic neutrinos is therefore the trigger for the onset 
of the cosmic acceleration. Such transition occurs at relatively late times (z m ~ 4 for 
the specific model considered here) thereby providing a possible solution to the cosmic 
coincidence problem. On the other hand, the oscillations of the DE scalar field around 
the minimum of the potential determine in turn oscillations of the average neutrino 
mass according to Eq. [T] that have very significant implications for the evolution of 
perturbations in the context of Growing Neutrino cosmologies, and that might provide 
a direct way to test, constrain, or disproof the model, as we will discuss in full detail in 
the next Sections. 



2. N-Body simulations 

The linear evolution of perturbations in the context of Growing Neutrino scenarios has 
been extensively studied by several authors (see e.g. |Mota et aT|2008t [Pettorino et al. 



2010| ), and some attempts have also been made to extend such investigations to the 



nonlinear regime by means of suitably modified spherical collapse numerical codes 
(e.g. by |Winterger st et aL]|2010| ). A common prediction of all these studies is that the 
neutrino fluctuations grow very rapidly after the transition redshift z nr due to the fifth- 
force mediated by the scalar field which is a factor 20 1 (i.e. around 5 x 10 3 times for our 
specific model) stronger than the standard gravitational interaction. This rapid growth 
leads to the formation of very large neutrino structures at scales beyond ~ 10 Mpc//j 
that grow nonlinear at z ~ 2. After this time the linear approximation no longer holds, 
and a full nonlinear treatment of the system is therefore required. The attempt to use 
modified spherical collapse algorithms, however, only mildly improves the situation 
by allowing to follow the evolution of nonlinear bound structures up to the onset of 
virialization, which happens at z ~ 1-5 for neutrino halos of size ~ 10 Mpc. This 
further limitation is mainly due to another peculiar effect of interacting DE models, i.e. 
the presence of an extra friction term for the coupled particles arising from momentum 
conservation. After the onset of virialization the system can therefore be studied only 
by means of suitable N-body simulations, that must include in their implementation all 
the specific features of interacting DE models, and in particular of Growing Neutrino 
scenarios. To this end, we have made use of the specific modification by [Baldi et al 



2010[ ) of the widely used parallel N-body code GADGET (Springel 2005 ) to run a series 



of large N-body simulations for the specific Growing Neutrino model described above. 
By using such numerical setup we have been able to extend the study of nonlinear 
neutrino structures down to z = 1, and to investigate the growth and the abundance of 
such structures in a realistic cosmological realization. The limit of z — 1 is due to the 
fact that neutrino velocities grow very rapidly into the relativistic regime due to their 
very large accelerations, and the Newtonian approximation on which N-body codes are 
based becomes not appropriate for such large velocities. 



Nonlinear evolution of Growing Neutrino cosmologies 
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Figure 1. Neutrino structures in a box of 120 Mpc/h. See Baldi et al. ( 201 1\ for 
further details. 



3. Results 



The results of our numerical simulations have shown for the first time how large neu- 
trino overdensities are seeded by the underlying Cold Dark Matter gravitational po- 
tential wells at the transition redshift z nr and how these bound structures subsequently 
grow nonlinearly and merge to form isolated neutrino lumps at scales beyond 10 Mpc. 

3.1. Neutrino large-scale structures 

The evolution of the neutrino distribution at different redshifts in a periodic cosmolog- 
ical box of 120 Mpc/h aside is shown in Fig. [I] The rapid growth and merging history 
of neutrino lumps between z - 1.8 and z - 1.3 is clearly displayed by comparing the 
left and the central panels. However, the right panel shows a significantly lower number 
of neutrino halos if compared to the situation at z = 1-3 and correspondingly a large 
fraction of free neutrinos that fill the whole simulation volume and do not belong to 
any bound structure. This peculiar evolution represents one of the main findings of the 
present work and can be explained in terms of the rapid oscillations of the neutrino 
mass due to the oscillations of the DE scalar field <p around the minimum of its effective 
potential: the change in sign of the time derivative of the neutrino mass m v determined 
by the oscillations implies a sort of "pulsation" of the gravitational potential associated 
to any given neutrino overdensity. In addition to this, the change in sign of the extra 
friction term also contributes to trap neutrino particles in bound structures whenever 
m v > while it favors the escape of neutrino particles from bound objects whenever 
ihy < (see Bal di et al.||2011| for more details). We have therefore shown for the 
first time that nonlinear neutrino lumps forming in the context of Growing Neutrino 
cosmologies are expected to "pulsate" following the oscillations of the DE scalar field. 

3.2. Neutrino halo mass function 

In order to quantify the expected observational signatures of nonlinear neutrino lumps it 
is necessary to determine, besides the specific time evolution of any individual neutrino 
halo, also a realistic statistical sample of such structures, in order to estimate their 
volume and line-of-sight densities. This is another issue that previous investigations 
could not properly address, and that has found a first assessment only with the present 
work. We have in fact identified neutrino halos in our simulations by means of a suitable 
modification of the widely used Friends-of-Friends algorithm and we have therefore 
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Figure 2. The neutrino halo mass function at different redshifts. 

computed a neutrino halo mass function (analogous to the standard CDM halo mass 
function) over the whole range of masses covered by our numerical simulations. The 
results of this procedure are shown in Fig. [2] where the abundance of neutrino halos as 
a function of mass is plotted for different redshifts. Also this plot clearly shows how 
the abundance of neutrino halos first increases in time up to z ~ 1 .4 and subsequently 
decays again, due to the scalar field oscillations and the consequent oscillations of the 
neutrino mass. 

4. Conclusions 

We have presented the results of the first N-body simulations of structure formation for 
Growing Neutrino cosmologies. Our results have shown for the first time how large 
neutrino lumps forming in the context of such cosmological models show a "pulsation" 
due to the rapid oscillations of the DE scalar field around the minimum of its effective 
potential. We have also determined, for the first time, a neutrino halo mass function 
that represents a fundamental quantity in order to put reliable observational constraints 
on this class of cosmological models. 
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